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COMPUTER PROGRAM FOR ANALYSIS OF COUPLED-CAVITY 


TRAVELING-WAVE TUBES 
by Denis J. Connolly and Thomas A. O'Malley 
Lewis Research Center 

SUMMARY 

A flexible, accurate, large- signal computer program has been developed for the 
design of coupled- cavity traveling-wave tubes. The program is written in FORTRAN IV 
for an IBM 360/67 time-sharing system. The beam is described by a disk model, and 
the slow-wave structure by a sequence of cavities, or cells. The computational ap- 
proach is arranged so that each cavity may have geometrical or electrical parameters 
different from those of its neighbors. This allows the program user to simulate a tube 
of almost arbitrary complexity. Input and output couplers, severs, complicated veloc- 
ity tapers, and other features peculiar to one or a few cavities may be modeled by a 
correct choice of input data. The beam-wave interaction is handled by a new approach 
in which the radio-frequency fields are expanded in solutions to the transverse mag- 
netic wave equation. All significant space harmonics are retained. The program was 
' used to perform a design study of the traveling-wave tube developed for the Communi- 
cations Technology Satellite. Good agreement was obtained between the predictions of 
the program and the measured performance of the flight tube. 


INTRODUCTION 

Numerous works describing large-signal computer programs for traveling-wave 
tubes have appeared in the past 20 years. A fairly complete list is given in reference 1. 
Some significant works not mentioned in reference 1 are references 2 and 3. 

The present report describes a program tailored to the analysis of coupled- cavity 
traveling-wave tubes on an interactive time- sharing system. Approximations not nec- 
essary to keep central-processing-unit (CPU) time within reasonable bounds are 
avoided. Extreme flexibility is provided in the variety of tube features that can be 
modeled, since each cavity in the slow-wave structure has individually entered input 



parameters. The program can handle lumped or distributed severs, input and output 
couplers, cavity match details provided for cavities near the end of a stack, voltage- 
jump velocity resynchronization with an arbitrary number of discrete steps (up to 1 per 
cavity), and velocity taper designs of almost arbitrary complexity. Backward waves 
can be handled by an iterative procedure similar to that described in reference 3. 

The program was developed to provide an analytical tool for the evaluation of 
coupled- cavity traveling-wave- tube efficiency enhancement schemes. It also may be 
used for routine traveling-wave-tube design. 

Test cases were constructed for comparison of computed results witii the asymp- 
totic Pierce theory. Computed results were also compared with experimental measure- 
ments on the 200-watt traveling-wave tube of the Communications Technology Satellite. 
In both comparisons, agreement was quite good. 


MODEL OF COUPLED- CAVITY TEA VELING- WAVE TUBE 

Our model of a coupled- cavity traveling-wave tube is illustrated by figure 1. The 
beam is treated as a series of rigid but mutually penetrable disks. A complete de- 
scription of the beam trajectory is obtained by following the disks contained in a single 
beam wavelength (ref. 4). The number of disks per beam wavelength is an adjustable 
input parameter but is usually set at 24. The diameter of a beam disk (2rj^ in fig. 1) 
is allowed to take on up to three different values during the course of a calculation. 

This option is provided to accommodate the fact that traveling-wave-tube beams tend to 
expand in the saturation region of the tube. 

Similarly, the charge and mass per disk are allowed to diminish by preselected 
amounts at two preselected positions within the tube. This allows the program user to 
roughly estimate the effects on tube gain and efficiency of various amounts of beam in- 
terception by the structure. 

The tube body is treated as a conducting tunnel (radius a) divided axially into a 

th 

series of discrete cavities, or cells. In the center of each cavity (length of k cavity 
is denoted by Lj^) is a gap of length 21 -^. Impressed across the k“^ gap is a complex 
voltage (Symbols are defined in appendix A. ) 

In the absence of abeam, FloqueUs theorem (ref. 5) provides a simple known rela^ 
tionship among the various Vj^. For a forward propagating wave in the passband of a 
uniform structure, 


''k+l ® '^k 


( 1 ) 


2 



where a and ^ are known from cold-test measurements. The power flowing along 
the structure is given by the following equation (from ref, 6): 



( 2 ) 


where Zj^, the interaction impedance, also is known from cold-test measurements. 

The voltage wave in the presence of the beam is calculated by addition of a com- 
plex induced voltage to the propagating wave voltage; that is, 


V = e"" ^ ^ V + A V 

''k+1 ® ''k ^''k+1 


(3) 


where is an induced voltage to be defined later. 

The body of a coupled- cavity traveling- wave tube is supposedly a periodic struc- 
ture. If this were the case, all cavities would have the same L, Z, aL, /3L, and Z. 
However, a real tube has input and output couplers, severs, and perhaps velocity 
tapers and cavity match details, all of which require individual, nonperiodic treatment 
in the model. Thus, the properties of individual cavities are separately specified as 
input parameters in the program so that a great number of tube-design variations such 
as complicated velocity tapers or voltage-jump configurations can be easily modeled. 
Equation (3) is easily modified to handle such a quasiperiodic structure by the use of 
local values of a. and /3: 


'^k+l “ "^k.k+l^k 


+ AV, 


k+1 


(4) 


where 


'^k,k+l 



-(0!L)j^ -i(PL)k 


(5) 


The factor yj ensures that power flow is conserved. 

Floquet's theorem (eq. (1)) strictly applies only to an infinite-length periodic struc- 
ture. We are actually applying it to a finite-length, almost-periodic structure whose 
properties may vary slowly with position. Equations (4) and (5) should thus be consid- 
ered as approximations, analogous to the phase integral method as used in ionospheric 
radio- wave propagation studies (ref. 7) . 
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MATHEMATICAL FORMULATION 


Beam Disk Dynamics 


The relativistic energy of a particle can be written as 


E = 



+ qVo 


Because energy is conserved through the electron gun, 



= mQC^ +qVo 


( 6 ) 


(7) 


which gives directly 



For a Vq of approximately 10 kilovolts, the relativistic correction is about 1.5 per- 
cent. For analysis of a device whose operation depends on near synchronism of the 
electron beam and the circuit wave, a 1. 5-percent error in velocity was thou^t to be 
too large to ignore. Thus, relativistic dynamics were used to analyze the beam mo- 
tion. The velocity change due to an applied force is given by 


dt 




( 9 ) 


Equation (9) is the one- dimensional form of equation (16-58) of reference 8. 
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Space-charge forces. - In a frame moving along with the beam, the space-charge 
forces can be determined by the solution to an electrostatics problem. The field of a 
point charge in a cylindrical tunnel is given by the following equation (from ref. 9): 



(2 - 6go)e 




p, s 


|z| 


cos s{(p - 


<Pq) 


sP^^s%, sPQ> 
^P, sps+l^^p, 


( 10 ) 


where r^ is the beam tunnel radius, and 


Cl for s = 0 

I 0 for s 7 ^ 0 




y ( 11 ) 


-y 


and the charge is located at (pq, (Pq, 0). Consider abeam disk of radius rj^ and 
thickness^ tj^, centered at (pq = 0; = 0), with uniform charge density cr. With the 

use of equation (11), it is a straightforward exercise to show that the potential in the 
tunnel due to the beam disk is given by 


^disk 



_2 

“P-" 


2 


r z > 


D 


1_ 


cosh Pp qZ for I z I < — 




P=1 


( 12 ) 


^In this analysis, the disk thickness remains an integer submultiple of the initial 
beam wavelength (Uq/I). A physically more reasonable treatment would allow the disk 
thickness to vary with its own local velocity. This would, however, add another dimen- 
sion to the table in which space-charge forces are stored and would substantially in- 
crease computing time. We believe the difference in results would be insignificant. 
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and, further, the mutual repulsion force between two identical disks whose centers are 
separated by distance Zq is 

OO 

^ 

I ■ 

^p, o’^a^ 

/ 

p=l 


z, SC 


sgn Zq 
,222 


r 


2 sinh^i^PL^e'^P-O'^o' 


I ^0 I “ 




1 _ e'^'P. . 3-Pp, O^D for 1 Zq I ^ tj3 


(13) 


where q is the total charge of a disk and is given by 

q = 

and sgn Zq is the sign of Zq. The axial force is unchanged by a Lorentz transforma- 
tion. However, axial lengths in equations (10) to (13) should be those in the moving 
frame. Thus, equation (13) can be used to calculate the space- charge forces directly 
in the laboratory frame with the following length correction; 



This analysis assumes that the velocity spread of the particles is small (nonrelativis- 
tic). The accuracy would, therefore, deteriorate near the output of a relativistic, very 
high efficiency tube. 

Radio-frequency electric-field forces. - If the complex voltage on a given gap is 
known, the electric field in the nei^boring region can be written in an e?q)ansion of the 
form (from ref. 10) 
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_ -mV V 

2 sinh(mi) / , L lo^n’^a* 


iwt-^ z 

I n 


( 16 ) 


where 


C„(m,/Jn) = 


2Z[mZ sinh(mZ)cos(y3j^Z) cosh(mi)sin(/3^Z^ 


(mi)^ + 





and j3qL is the lowest order phase shift per cavity. The parameter m is chosen to 
best suit the geometry of the gap. In the case of very blunt beam tunnel tips, a zero 
value for m is appropriate, and equation (16) reduces to 


E 


z 



sin(/3^Z) Io(Tj^i’) -i/3j^z+io)t 

B I If.(y r ) 

'^n O' n a' 


(17) 


The expansion of equations (16) and (17) is appropriate to a forward- traveling circuit 
wave. An arbitrary combination of forward- and backward- traveling waves can be 
represented by 


oo 


E 


z 



n=-oo 


+ Vge 



iojt 

e 


(18) 


where Vj, and Vg are the complex amplitudes of the forward and backward voltage 
wave, respectively. For convenience, we have defined 
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(19) 


E^(r) = 


Ip(y^r) 

2L siiih(m?) Io('Viir^) 


Consider a beam disk of radius r^^, thickness tjj, and charge q, and centered at 
(r = 0; z = Zq). a straightforward calculation shows that the force on the disk due to 
the radio-frequency (RF) electric field is given by 


/3„t 


qm 


L sinh(mZ) 


Cn(m,^n)^l<ynV 


sin • 


n'D 


Pn- 
^ 2 


/^e 


z. 


lVj,e 


n=-oo 


( 20 ) 

For use in the program, many of the parameters in equation (20) take on a subscript k, 
since they may vary from one cavity to another. Thus, 


qm 


z,k, RF sinh(mij^) 


'^n.k<”-l^n.k>^lO'n.kV 

■>'n,k'’b'o<'''n,k''a> 




’n,k' 


"D 


n=-oo 


X /^e 


V Lict 

Vf k® ^'^B.k® /® 


( 21 ) 


where Zj^ is measured from the center of the k*^ g^, and tj^ is the disk thickness 
in the laboratory frame. 

Voltage-jump electric- field forces . - If a dc voltage Vj is imposed across a gap 
in the slow-wave structure, the resultant electric field can be analyzed in a manner 
similar to the analysis of the radio- frequency electric field discussed earlier. We first 
note that the axial electric field is an even function of z. We also note that the geom- 
etry of the problem typically satisfies the criteria L » r^, and L » 1. Hence, the 
fields near the gap have decayed to a fairly low level by z = ±(L/2) . We may therefore 
expect that little error will be introduced if we expand the dc gap field in a Fourier sine 


8 



series over the finite interval -(L/2) < z < (L/2). The potential function which satis- 
fies LaPlace's equation and the boundary condition at ±(L/2) and has the required be- 
havior at r = 0 is 


(p = Aq + B(jZ + ^ sin(k^z)lQ(k^r) 
n=l 


with 


k 

- L 

For our purposes, Aq is arbitrary, and since 



we have 



n=l 


sin(k^z)lQ(k^r) 


and 


E 


z 



B k 
n n 


cos(k^z)lQ(k^r) 


The coefficients B^^^ are available to fit the assumed field shape in the gap which 
take to be of the form 


E^(r-r^) 




n=l 


r 


Eq cosh mz if \z\ I 


if Izl > i 


■N 




( 22 ) 


(23) 


(24) 


(25) 

we 


(26) 
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where 


The boundary condition at 


mV, 


Eq = 


2 sinh(mZ) 


r = r leads in the usual way to 

3. 


B k = 
n n 


mVjC^(m,kn) 

L sinh(mZ)lQ(k^r^) 


(27) 


(28) 


with 


C„(m.k„)=2 


m sinh(mZ) cos(k^i) + cosh(mZ) sin^^i) 

2^,2 

m i-k„ 


( 29 ) 


Thus, equation (25) becomes 


mV, 


E =--i- 


L L sinh(mi) 


oo 


n=l 


(30) 


A strai^tforward analysis shows that the force exerted by this field on a beam disk is 
given by 


r 



+ 2m 

QQ 

\ <^n V Vo> — 

L 

sinh(mZ) 

i 




n=l 


If IZqI s- 


L - t,. 



f 

oo 

sin 


L 

2tp sinh(m^) j 

n=l 

Vd 

- 2 J 

-> 


Obviously, equals zero if |Zq| > (L+tj^/2. 


L - t. 


D . 


IZnI 


2 


( 31 ) 
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Beam-Wave Coupling Analysis 

Let E^, and Eg, Hg be solutions to Maxwell’s equations for current sources 
3^ and Jg, respectively. It is easily shown that 




d«e 


(32) 


where n is an outward-directed unit normal, and S is a closed surface ho unding the 
volume and asterisk denotes complex conjugate. We shall use equation (32) to de- 
rive an expression for the fields radiated by a current source in a coupled- cavity slow- 
wave structure. Let ** denote a section (volume) of a slow- wave structure bounded by 
the outer, cylindrical, conducting wall Sq and any two planes Sg and Sj, normal to 
the axis (see fig. 2). It is immediately clear that the integration over Sq makes no 
contribution to the left side of equation (32), since the tangential electric field must 
vanish on a conducting surface. Thus, equation (32) can be written as 



= - ^ [(e^ • h) + (®2 • »3) 

Orthogonality relation for periodic structures. - In equation (33), let Ej^, be a 
backward-propagating wave Eq g, Hq g whose source is outside the volume of inter- 
est. Let Eg, Hg be a forward-propagating wave Eq g, Hq ^ whose source is also 
outside the volume of interest. The volume integral for this case is zero, and equa- 
tion (33) can be written as 

[(®0, B ^ ^0, f) (®0, F ^ ^0, b)] ' ®z " y [(®0, B ^ ^0, f) 

Sf Sg 

+ (\f^Vb)] • = K (34) 


where K is some constant independent of the axial position at which the surface inte- 
gral is taken. 


11 



are 


To determine the value of K, consider the special case in which Sg and Sp 
separated by exactly one periodic length. Since the region is source free, the fields 
must obey Floquet’s theorem and hence 





Excepting the trivial case where jSL is a multiple of tt, equations (37) and (34) re- 
quire that K = 0. 

Fields radiated by a current source in a slow-wave structure . - Let a source 

be present in a sin^e cavity of a periodic coupled-cavity chain. The 
surfaces Sg and Sp are the axial boundaries of the cavity (z = ± L/2). Let 
Er> T 5 , Hn Ti be the fields of a backward wave whose source is outside the cavity ** and 
whose amplitude is chosen to represent unit power flow. Similarly, let Eq p, Hq p 
be a unit- amplitude source-free forward wave. Outside the cavity ««, the fields ra- 
diated by J„ within if are given by 

a. 





r ^ 
^F ®0, F 


if z > — 


^F ^0,F, 


(38) 
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( 39 ) 


a ^®0,B 



.. L 

if z < 




where a^, and are complex constants to be determined. 

In equation (32), let Ej^, be replaced by Eq g, g, and let Eg, Hg, J£ 


be replaced by E„, H^, J^. The result is 

a. a. a. 


; /* [(®0, B ^ ^0, b) (®0, B ^ ^0, b)] ■ S f [(®0, B ^ f) 


+ (Eq.F ^ ^0,b)1 ■ “ / ^0, 


g ■ d» (40) 


The integral over Sg equals 4, because Eq g, Hq g have been normalized to unit 
power flow. The integral over Sp equals 0 by virtue of equation (34). Hence, 


a. 


B 


-;A 


b' 


(41) 


A similar analysis with E^, in equation (32) replaced by Eq j,, Hq g leads to 


^F " 


f 


J„ d«* 
a 


(42) 


The fields Eq g, Eq g are expressed in the form given by equation (18), The gap 
voltages Vq g and Vq g must satisfy (|V|^/2Z) = 1, since we have specified unity 
power flow in the source- free fields. For convenience, we specify zero phase for 
Vq g and Vq g, so that 


and, therefore. 
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( 44 ) 


and 



(45) 


Representation of radiate_d fields by induced voltage. - Equations (38), (39), (41), 
and (42) describe the fields radiated by J„ outside the cavity that contains J . The 
result may conveniently be represented by an induced gap voltage AVj, which propa- 
gates in the forward direction and another, AVg, which propagates in the backward di- 
rection. With phase referred to the gap of the cavity containing J , the induced volt> 

a 

ages are given by 


and 


AV. 


B 



(46) 


(47) 


Equations (41) and (42) can be used to calculate the effect of J on the circuit wave 

cL 

everywhere except within the cavity that contains J„. These results are not valid 

a, 

within the cavity, because equations (38) and (39) do not ^ply there. To obtain a result 
valid within the cavity, consider figure 3, which depicts the power balance within the 
cavity. The power, P„, added by J is, evidently, 

3, a. 


14 



p =— 
® 2Z 


\Re (\ 


2Ke (VpAVp + V* AVg 


+lAVj,|^ + |AVgl^ 


■^Re 

z 



AVj, + V3 + 



AV. 


(48) 


Or, with the use of equations (46) and (47) , 


P_ =-A 


/ \ 

/ \* • / \* 


/ AV \ 1/3 z / AV„\ -i^ z 

. /j e ” * , nj 


d«* (49) 


But, we also know that P is given by 

ci 


Pa = 


Re / 

2 Jn 


E • J dt* 
a 


(50) 


Hence, correct results will be obtained if the forward- and backward- propagating 

gap voltages which act upon the beam within the cavity are set equal to jvj, + (AVj,/2)J 

and (AVg/2^, respectively. In the preceding analysis, Vj, and are the gap 

voltages due to fields which propagate in from neighboring cavities and are independent 
of J^; the induced voltages AVp and AVg are given by equations (46) and (47) . 

Power-Balance Verification 


Consider the section of a traveling-wave tube beginning at the input end of the tube 
and terminating at the exit plane of the k cavity. At the input end, power enters in 
the form of beam kinetic energy (IqVq) and input radio-frequency drive (Pq) and leaves 
in the form of backward wave power (Pg j^) propagating back into the driving source. 

At the output end of the test section under consideration, power enters in the form of 
backward-wave power (Pg propagating back from cavity k+1 and leaves in the 
form of beam kinetic power (Pj^ j^) and forward-wave radio-frequency power (Pg ^^). If 
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"there are any voltage-jump cavities within the test section, they add power 

k 

j=l 

Finally, power is dissipated by joule heating within the test section 

k 

2^^LOSS,j 

j=l 


A test parameter can be defined by writing the power balance within the test section. 
In the notation defined above. 


k k 

^0^0 ^0 " ^B, 1 k+1 " ^b, k " ^F, k ^J, j " ^ ^LOSS, j 

j=l j=l 


or 


a 


PB,k 


^b,k ^F,k ^B, 1 ^LOSS, j " ^J, j “ ^0 " ^B,k+1 

( 52 ) 


IqVo 


The parameter “pg is calculated for every cavity and printed out as an independent 
check on other calculations; obviously, it should equal unity. Individual contributions 
to this parameter are 




( 53 ) 
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where is the rest mass of a beam disk, u. is the axial velocity of the i*"^ disk at 
the exit plane of the k cavily, and the sum is over the disks in an entering beam 
wavelength. Also, 


^LOSS,j 


^b, 0 ^0^0 


p 

22 ^ 


IV, 


B,k 


B,k' 


2Z, 


Pj, j ' 




1 - 




(54) 

(55) 

(56) 

(57) 

(58) 


Equation (58) is correct only if the forward and backward waves are 90*^ out of phase 
and thus not coupled by cavity wall losses. However, the errors due to use of equa- 
tion (58) will tend to average out to zero when a large number of cavities with a variety 
of phase angles are included in the test section. Usually, P-p . » Pp ., which further 
diminishes the importance of this error. Finally, the joule-heating loss is usually not 
one of the largest terms in equation (51) . Therefore, equation (58) seems to be a fairly 
safe approximation. 

A more important error is neglect of space- charge potential energy transported by 
the beam. For hi^ perveance devices, this may equal a significant fraction of the 
beam kinetic power (ref. 11). However, the magnitude and direction of the error can 

— q /O 

be estimated. In the case of low perveance devices (« O.lxlO"” AV^'n, it is negli- 
gible. And the energy balance parameter serves only as a running check on the other 
calculations. Because accounting for beam potential energy would significantly in- 
crease computer time, this was not done. 
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PROGRAM STRUCTURE 


Input Data 


The Inputs required by the program, the units in which the inputs are entered into 
the program, and the symbol used for each input are as follows: 


Beam voltage, V 

Beam current, mA 

Input power, dBm 

Frequency, GHz 

Beam tunnel radius, cm 

r beam radius, cm 

Electric field shaping parameter 

Number of beam radii 

Number of cavities 

Number of disks per beam wavelength 

Number of nodes per cavity 

For the cavity: 

Cavity periodic length, cm 

Gap length, cm 

Propagation loss across cavity in forward direction, dB . 
Propagation loss across cavity in backward direction, dB 

Interaction impedance, ohms 

Phase shift parameter 

dc voltage jump, V 


Vq 

.... Iq 

. . . Pq 
. . . . f 

' ■ ‘ ^a 

• • %,i 

... m 

. . . N^ 
. . . N^ 
. . . N^ 

• • • N, 

• • • 

. . 21 ^ 

. . . . Zk 

[(gjL)/,]^ 

• • • ''j.k 


Beam Disk Dynamics 

The computational problem discussed in this section is how to calculate beam disk 
trajectories in a suitable manner. Ordinarily, the computational procedure is to solve 
for the disk acceleration, as given by equation (9). Numerical integration of this equa- 
ls 


I 



tion then yields u. (t) , the velocity of the i disk at time t. Numerical integration of 
u.(t) yields z.(t), the axial position of the i^ disk at time t. Instead of using the 
functions u.(t) and z.(t) to describe the trajectories, it is more convenient for our pur- 
poses to let z be the independent variable and to calculate the functions dtj(z)/ dz and 
t-(z). Here, dt.(z)/dz is the reciprocal of the velocity of the i^ disk when it reaches 
the axial position z, and t.(z) is the time of arrival of the i disk at the axial posi- 
tion z. In terms of these new functions, the equations of motion are 


3/2 



dz 


2 



m 


D 




for i= 1, 2, 



(59) 


til 

where F_ . is the total axial force on the i disk. 

Z y 1 

The coordinate system used in the program has its origin at the beginning of the 
first cavity. The k*^“ cavity is divided into equal parts of length AZj^ given by 


Az, = 


(60) 


.th 


The n node in the first cavity, denoted by z^ is defined by 


z^ ^=aAz^ for n = 0, 1, 


( 61 ) 


The n*'^ node in the k^^ cavity, denoted by is defined by 


k-1 


^n,k 




for" 


m=l 


n = 0, 1, ..., N 


k = 2, 3, ..., N 


(62) 


A point which lies midway between two nodes z^ and is denoted by 

Such a point will be called the n^ antinode in cavity k. The numerical integration 

2 2 

of equation (59) is such that the functions d t^/dz are evaluated at nodes, and the 
functions dt^/dz and t.(z) are evaluated at antinodes. 

At some point in cavity k, we have trajectory data at node n-1 and antinode n-1; 
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2 2 

that is, the functions d t.(z„ , ,)/dz , dt.(z„ i dz, and t.(z„ i ,) are known. 

1 X X IX“*2 > 

2 2 

The problem now is to obtain d tj(z^ j^)/dz , dt.(z^^ ^'i^^n+l k^* 

first obtain approximations for dt.(z^ ^i^^n k^ equations 


( 63 ) 


dz 


ti(z„,k> = * 


.,z . ^^.,1 A<vi.k> /^\\ 

^“"2,k 2 2 .2 \ 2 / 

dz ' ' 


(64) 


Knowing, at least approximately, dt.(z O/dz and t.(z ,), we can evaluate 
2 2 

d t.(z .)/dz from equation (59). We now obtain dt.(z„.i ,j/dz and t.(z^,i . ) by 

X II, XV X II**" 2 5 ^ ^ H“*" 2 5”^ 

the equations 




dz 


dz 




dz 


(65) 


k^ ‘^^^i^^n k^ ^4 

t-(z 1 ) = t.(z 1 1 ) +■■■ Az, + LJli^ ± 

i' n+4,k' i' n-t,k' “Tc 


dz 


dz 


( 66 ) 


Space- charge forces . - The space- charge force on one disk due to another with a 
separation distance of Zq is given by equation (13). Use of this equation for calcu- 
lating space- charge forces on each integration step would be intolerably slow. The 

method used here is to set up N, tables, where each table contains 50 data points 

o th 

corresponding to 50 separation distances. The n data point would be the space- 
charge force between two disks, as calculated by equation (13), separated by the n*"^ 
separation distance. In the computer simulation, the space-charge force on one disk 
due to another is calculated by the use of linear interpolation on the appropriate table. 

Fortunately, the space-charge force on a given disk is a slowly-varying function of 
the independent variable z. This is especially true in the beginning cavities of the 
tube. To take advantage of the slowly- varying nature of the space-charge force, the 
calculation is not done on every integration step. The calculation of space-charge fore- 
es is initially done eveiy step, where M is an input. During the simulation. 
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M may change in a manner to be described later. On those integration steps where 

space-charge forces are not calculated, they are approximated by quadratic curve-fits. 

Each disk has its own curve-fit. These curve- fits are obtained by fitting a quadratic 

polynomial through the latest three calculations of space-charge forces. Let us say 

that these three calculations were done when z was equal to z„, z. , and z^. At an 

a’ b’ c 

integration step where a new space-charge force calculation is due, at z = z^, the val- 
ues calculated are compared with those obtained when the curve-fits based on the points 
z^, Zj^, and z^ are evaluated at z^. If the difference in the values for any disk is ex- 
ceeded by a specified tolerance, the next calculation of space-charge forces occurs 
steps later. In general, M is halved each time the specified tolerance is exceeded. 

For calculating space-charge forces, the positions of all disks must be known at a 
time t. This information is not readily available, since z rather than t is the inde- 
pendent variable. To calculate the space-charge force on the i*^“ disk located at some 
z, we must determine the positions of all the other disks at the time t.(z). Consider 

the j disk. The time of arrival, t.(z), at z is known. By the periodicity of the mo- 

J fW 

tion, the disk which originally was m wavelengths behind the j disk arrives at z 

at the time t. (z)+mT. Similarly, the disk which originally was m wavelengths ahead 
th ^ 

of the j disk arrives at z at the time L(z) + (-m)T. We consider all such disks 
and find m such that 


j t.(z) - [tj(z) + m"^ j = minimum (67) 

th til 

We may call the disk in question the j disk in the m*^ cycle. We assume that the 
disk in any other cycle is too far away from the i*'^ disk at the time t.(z) to make 
a significant contribution to the space-charge force on the disk. Let z. be the 
position of the disk in the m*"^ cycle at the time t.(z). We obtain z. by assuming 
this disk has constant velocity Uj(z) in the time interval spanned by t.(z) and 
tj(z)+mT. Then z. is given by 


Zj. = z + m(z) |t.(z) - [tj(z) + (68) 

The space-charge force is then calculated on the basis of the separation distance 

z-z.. The calculation is done for all disks, excluding the ,i*^ disk, and the results 
3 th 

are summed to obtain the total space-charge force on the i disk. When the separa- 
tion distance ^ - Zj is greater than a specified number, the space-charge calculation 
is not done for this term in the summation. 

Rad io- frequency electric-field forces . - The radio- frequency electric-field force 
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on a disk is given by equation (21) . The in this equation is the axial position rela- 
tive to the center of the cavity. In the present notation, z, wHl be replaced by 

fVi fV* 

y , . Thus, y is the axial position of the n*^“ node in the k“ cavity relative to 

iij K Hj k 

the center of the cavity. In equation form, y^^ is related to z^ by 


We define by 


'k-1 


y i_ = z t ” I \ L I 

•^n,k n, k I / m g 


m=l 


(69) 
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’m,k''D 


m=-oo 


1^1 0, 1, . . • , 

for<^ (70) 

[k=l, 2, ..., 


fVi 

where jjL = Then, the electric field force on the i*^“ disk located at is 

given by 






k^F, k ^n, k^B 


,0' 




(71) 


The Qjj are calculated and stored in a table when cavity k is entered. The number 
of terms of the summation used to calculate each table entry is an adjustable input pa- 
rameter. If cavity k has the same Lj^, /3 q as cavity k-1, then 
and a new table does not have to be calculated. When a new table does have to be cal- 
culated, the old table is no longer needed, and the new table may occupy the same stor- 
age space as the old table had occupied. The electric- field force is evaluated by equa- 
tion (71). 
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Voltage- jump electric- field forces . - The voltage-jump electric- field force on a 
disk is given by equation (31) . The technique for evaluating voltage-jump electric-field 
forces is essentially the same as for evaluating radio-frequency electric-field forces. 
We define R . by the equations 




,jVd' 


1 +- 


hL 




Ij^ sinh(/j) 


m=l 









sinh(ji) 



k tj., 
m JJ 








m=l 


if < ^k^^D (72b) 

2 2 


Let Vj be the voltage jump associated with cavity k. Then, the electric-field force 
on a disk located at z , is given by 


^n, k ^n, k^J, k 


(73) 


The are calculated and stored in a table and have the same features as the 

for radio-frequency electric-field forces. ’ 
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Beam-Wave Ck>upling Analysis 


The parameters of interest in the beam-wave interaction are the induced voltages, 
AVp and AVg. These are obtained from equations (46) and (47) . We proceed to use 
these equations to develop expressions for AVj, and AVg which are computationally 
convenient. The appropriate Fourier component of the beam current, Jj^(z), is given by 


•'t-T 




We define by 




With the use of equations (19), (46), and (47), we obtain for cavity k 



(74) 


(75) 


(76) 


AV. 


B,k 




m,k^ 


Jj^(z)dz 


(77) 


In equations (76) and (77), the integration variable z is the axial position relative to 
the center of cavity k. The integrals in equations (76) and (77) are evaluated with the 
use of Simpson’s rule. Thus, 


N_ 


z OO 

n=0 m=-«> 


^m,k® 


^^m,k^n,k 




(78) 
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where 


N. 




Z-.Z 

n=0 m=-oo 


P , J (z ) 

^m,k® 


— ifn = Oorn=N 


6 s J— ifnj^O, N , and n is an odd integer 


n 


— if n 0, n^N , and n is an even integer 


We define by 


( 79 ) 






m,k^ 
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, k^n, k 


( 80 ) 


The values are calculated and stored in a tsble that has the same features as the 

Qjj and tables. Appendix B gives a derivation of the following expression for 




( 81 ) 


Our final ej^ressions for AVp and AVg ^ are 
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(82) 


N„ 


^''F.k = 2 «n®n,k‘'l<Vk> 


n=0 


N 




n=0 


(83) 


Calculating the effects of the backward wave requires an iterative procedure. In 
the first pass through the tube, the AV^ values are calculated and stored for each 
cavity k. From the values of AVg tiie backward voltages for a cavity chain from 
cavity kl to cavity k2 are obtained from equations (4) and (5) and the results of the 
Beam-Wave Coupling Analysis section. 


= -AV^ 


^B,k2 ““’'B,k2 


^B,k — ^ rB,k+i ^ — 



\ -(o;L)b k ■^<3iB)k 

k+l 


(84) 


for k = (k2 - 1), (k2 - 2), , kl (85) 


With the known values of Vg for the desired cavities, a second pass is made throu^ 
the tube. The second pass yields a new set of AVg values, which can be used for 
calculating a set of Vg ^ values for a third pass. The process continues until conver- 
gence is obtained. In many applications, the iterative procedure will converge faster 
if the iteration is done on one tube section at a time. This would be true, for example, 
for a tube having severs. The program has the capability of performing the iterative 
procedure for an arbitrary set of cavity chains. 

The computational procedure in the beam-wave interaction process is as follows: 

(1) When cavity k is entered, the summation equation (80) is calculated 

and stored for each node n. If cavity k has the same Lj^, and r^ as has 

cavity k-1, then S = S . , , and a new table does not have to be calculated. 

(2) A first ^proximation to AVg is obtained by assuming the disks have con- 
stant velocity in cavity k. Thus, 
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<Vk- ^0,k> 


( 86 ) 


dti(Zo 


With t.(z^ j^) known, j^) can be calculated from equation (81) , and then AVj, 

can be calculated from equation (82) . 

(3) Let the forward voltage for cavity k be 
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(87) 


If there is a backward voltage, it is assumed to be known from use of equation (85) with 
the set of AVg obtained from the previous pass through the tube. 

(4) Numeric^ integration of the equations of motion is used to obtain j^) for 

n = 1 to n = N„ . 

fVi 

(5) At the n “ node in the cavity, after tj(z^^ j^) has been obtained, k^^l^^n k^ 

jjj ^ ^ 

and k'^l^^n k^ calculated and stored. These quantities are the n “ terms in 

the summations of equations (82) and (83) . 

(6) When the last integration step in cavity k has been done, abetter approxima- 
tion of AVp is calculated from equation (82). The AVp so obtained replaces the 
old AVj, calculated in step 2. Also, AVg is calculated from equation (83). The 
AVg values are stored for use in the next pass through the tube. 

(7) If additional accuracy is required, steps 3 to 6 are repeated. 

(8) Steps 1 to 7 are repeated for cavity k+1. 

(9) When the last cavity is done, and the effects of the backward wave are to be de- 
termined, a second pass is made through the tube. The backward voltages for the sec- 
ond pass are determined from the AVg values of the first pass. 

(10) As many passes are made through the tube as are required to obtain conver- 
gence. 


Power-Balance Verification 

The equations for power balance, equations (52) to (58), are in a form that can be 
implemented without any modifications. It is v/orth noting that one of the advantages of 
using z rather than t as the independent vari^le is in calculating the kinetic power, 
P,, From equation (53), P. , is easily calculated, since the values of u.(z) are 
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known. If t were the independent variable, one would have to keep track of each disk 
as it passed the point z in question and calculate u.(z) by an interpolation procedure. 


Summary of Program Steps 

The entire computational procedure may be summarized by the following steps: 

(1) Input data are read. 

(2) The tables for spaco- charge forces are calculated. 

(3) Numerical integration of equations of motion is begun with z as the indepen- 
dent variable. At each step, the radio- frequency forces and, if required, the voltage- 
jump forces are calculated. The space-charge forces are calculated either by quadrat- 
ic curve fits or by linear interpolations on the appropriate tables. 

(4) At the end of each step, it is determined whether data are to be stored for 
plotting. 

(5) When a new cavity is entered, let us say the k*^ cavity, the data for cavity 

k-1 are printed out. If the parameters of cavity k are different from those of cavity 
k-1, the tables which are required for evaluation of radio-frequency forces, voltage- 
jump forces, and induced voltages are calculated. A first approximation of AV^, is 
obtained by assuming the disks have constant velocity throughout cavity k. The forward 
voltage for cavity k is obtained by attenuating and phase- shifting Vj, and adding 

the result vectorially to AVp j^2. 

(6) When the last integration step in cavity k has been done, a better approxima- 
tion of AVp is calculated. Also, AVg is calculated. 

(7) If additional accuracy is required, a second pass is made throu^ cavity k. 

(8) For cavity k+1, steps 5 to 7 are repeated. 

(9) If the effects of the backward wave are to be determined, a second pass is made 
through the tube. 

(10) As many passes are made through the tube as are required to obtain conver- 
gence. 


RESULTS AND DISCUSSION 

The coupled- cavity program has been tested by appropriate comparison with the 
asymptotic Pierce theory (ref. 12) and also with experimental results obtained on the 
200-watt traveling- wave tube of the Communications Technology Satellite (refs. 13 
to 15). 
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Comparison of Analytical Results With Asymptotic Pierce Theory 


In program runs to be compared with the asymptotic Pierce theory, the relativistic 
corrections and all space harmonics other than the first forward one were suppressed. 
The backward wave is, of course, ignored in both calculations. Program inputs de- 
scribing circuit and beam parameters were chosen to achieve predetermined values of 
Pierce’s C and QC over a range of values of b. Input power and tube length were 
chosen so that, for each value of C and QC, gain exceeded 20 decibels over most or 
all of the range of b, and tube output remained well below saturation. In all runs, the 
step-size parameters were chosen as 24 disks per beam wavelength, and 16 axial mesh 
points per period. The space- charge calculation was cut off at a disk separation dis- 
tance of 6. 5\g/24. 

In each case, tube gain, calculated as above, was compared with values given by 
the asymptotic Pierce theory (i. e. , from ref. 12). 

G = A + BCN (88) 


A = 20 log 
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N= (91) 
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B = 54.6 • (6j^) (92) 


QC = 



2 .i.T> 2 

Cl) + Rc»)p 


(93) 


The space-charge reduction parameter R was obtained from figure 9. 3-3 of refer- 
ence 12 which was stored as a table. 
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At first, the propagation constants 6 ^^, Sg. and 63 were obtained by calculating 
the roots of the cubic determinantal equation appropriate to the small C approxima- 
tion (ref. 12 ). 


5 ^= 4QC (94) 

(-b + id + i5) 

However, at moderate values of ICb| (i.e., when lCb| >0.1), rather large discrep- 
ancies appeared between the results of the asymptotic Pierce theory and the coupled- 
cavity program. The discrepancies did not significantly diminish when the number of 
disks per beam wavelength or mesh points per period was increased. It appeared, 
therefore, that the fault might be in the approximations leading to equation (94). A 
more exact fourth^order determinantal equation was then obtained from equation 
(10.1-35) in reference 12. 

9 (1 + Cb - iCd)(l + Cb)^ 

6 = — 4 QC (95) 

(-b + i6 + id) - — - (b - id)^ 

2 L J 


It is clear that equation (95) reduces to equation (94) if one assumes C « 1 . However, 
2 3 

the error in 6 will be of order (1 + Cb) - 1 « 3Cb. Subsequent calculations of 
asymptotic Pierce gain used the three roots of equation (95) corresponding to forward 
waves. Roots of equation (95) were found by a simple iterative procedure that used a 
root of equation (94) as a first approximation. 

After this refinement, agreement between the coupled- cavity program and the 
Pierce asymptotic theory was found to be satisfactory. The comparison for a range of 
parameters is shown in figure 4. 


Comparison of Analytical Results With Measured Data 

Tube design . - The design of the 200-watt traveling-wave tube for the Communica- 
tions Technology Satellite (CTS) was determined largely by others, prior to the develop- 
ment of our large- sigpal program. Subsequently, because of project schedule restric- 
tions, only a few minor design changes which addressed major defects were allowed. 

The most important of these changes are described in reference 15. 

The design beam voltage, beam current, nominal beam diameter, and tunnel diam- 
eter are 11.1 kilovolts, 70 milliamperes, 0.635 millimeter, and 1.27 millimeters. 
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respectively. The standard cavities in the gain sections of the tube have a periodic 
length of 3.175 millimeters and a gap length of 0. 965 millimeter. 

The slow-wave structure has a rather large cold bandwidth of about 1300 mega- 
hertz, with the 8 5- megahertz CTS band centered near a phase shift per cavity of 1. 2 tt 
in all sections of the tube. 

The number and t57pes of cavities in the CTS tube are shown in figure 5. The three 
gain sections are separated by severs. The three velocity- taper sections are separated 
by transition regions in which the phase velocity (or periodic length) changes gradually 
from one cavity to the next. Each sever consists of two modified cavities, each con- 
taining tuning stubs and silicon carbide - beryllium oxide loads. Forward and backward 
circuit waves are completely terminated at a sever, since there is no coupling slot in the 
wall between the two sever cavities. Two additional cavities adjacent to the input and 
output couplers and on each side of the two severs are modified for match purposes. 
These cavities have extra long coupling slots and thus extra wide cold bandwidth; they 
provide a reduced contribution to the gain of the tube and have much different electrical 
properties than their nei^bors. 

The beam tunnel tips are rather blunt, so the parameter m in equation (16) was 
taken as zero. 

Raw cold-test data obtained with uniform stacks of 5 to 10 cavities were available 
for many of the particular cavity geometries including all of those for which there was 
more than one cavity per tube. Phase velocity and interaction impedance, at the fre- 
quencies of interest, were calculated from the cold-test data by use of published meth- 
ods (refs. 6 and 16). Results were interpolated or, in two cases, extrapolated to cover 
those lO cavities for which cold- test data were not available. Accurate insertion loss 
measurements on brazed tube sections were also available. Values for propagation 
loss per cavity were obtained by dividing the insertion loss measurements, in a reason- 
able way, among the cavities within a given section. Resulting input data used for the 
calculation are given in tables I to IV. 

Comparison of results . - Computed results and corresponding experimental data 
for the CTS traveling-wave tube are presented in figures 6 to 8. All measurements 
were obtained with abeam voltage of 11.3 kilovolts and beam current of 71 milliam- 
peres. Beam current interception was t5q}ically 1 milliampere under small- signal con- 
ditions and 4. 5 milliamperes at saturation conditions. Calculations were also made 
with a beam voltage of 11.3 kilovolts, but the beam current was taken as 70 milliam- 
peres in order to make an approximate allowance for current and power removed from 
the beam by interception. The beam diameter used in the calculation was tiie nominal 
value of 0.635 millimeter (one-half the tunnel diameter). Results are shown for the 
first iteration with backward wave ignored, and for the final result after several itera^ 
tions with the backward wave included in the calculation. Three iterations were usually 
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sufficient to converge to the point where further change was negligible. 

Figure 6 shows small- signal gain as a function of frequency. Ignoring the gain 
ripple, which is assumed to be due to internal reflections, the agreement between mea- 
surement and calculation is excellent. The first iteration, without considering the 
backward wave, is also in good agreement with experiment. However, a significant 
discrepancy begins to appear near the lower edge of the band. In this particular tube, 
the IT point was around 11. 925 gigahertz. 

Figure 7 shows saturation drive as a function of frequency. Again, agreement be- 
tween calculation and experiment is excellent. The first iteration again shows appre- 
ciable errors only near the lower edge of the frequency band. 

Figure 8 shows saturation output power as a function of frequency. Agreement be- 
tween calculation and experiment is within about 0. 5 decibel (about 3 percent of input 
beam power), with the calculated power usually being somewhat hi^er than the mea- 
sured power. This is comparable to the amount of additional power dissipated by beam 
interception under saturation conditions which also varied somewhat with frequency. 

The effect of ignoring the backward wave was a further increase of 0. 5 to 1 decibel in 
computed output power. This discrepancy was observed everywhere in the passband 
and became a bit larger near the lower edge of the band. 

Limited calculations were also made with other values of nominal beam diameter. 

A 10-percent increase in beam diameter produced about a 1-decibel increase in gain 
and a 0. 1-decibel increase in output power; a like decrease in beam diameter had the 
opposite effect. No appreciable change was observed in the shapes of the various 
curves as beam diameter was varied. In all of these calculations, the internal check on 
power balance was satisfied extremely well, with the discr^ancy rarely exceeding 
0. 1 percent and never exceeding 0.2 percent. Computer central processing unit (CPU) 
time required for a sin^e iteration was less than 1. 5 seconds per cavity. 


CONCLUDING REMARKS 

A flexible, accurate, large- signal computer program has been developed for analy- 
sis of coupled- cavity traveling-wave tubes. The program was written in FORTRAN IV 
for an IBM 360/67 time- sharing system. The central processing unit (CPU) time for a 
single pass through the tube is about 1. 5 seconds per cavity. The beam is described by 
a disk model, and the slow-wave structure by a sequence of cavities, or cells. The 
computational approach is arranged so that each cavity may have geometrical or elec- 
trical parameters different from those of its neighbors. This allows the program user 
to simulate a tube of almost arbitrary complexity. Input and output couplers, severs, 
complicated velocity tapers, and other features peculiar to one or a few cavities may be 
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modeled by a correct choice of input parameters. The beam-wave interaction is 
handled by a new approach in which the radio-frequency fields are expanded in solutions 
to the transverse magnetic wave equation retaining all significant space harmonics. 

The program was used to perform a design study of the traveling- wave tube developed 
for the Communications Technology Satellite (CTS) . The computed results were in 
good agreement with the measured data from the flight tube. The agreement for satu- 
ration drive and small- signal gain was within 1 or 2 decibels, and saturation output 
power agreed within about 3 percent of input beam power. The internal check on power 
balance was satisfied within 0.2 percent . 

The small- signal gain calculated by the program was also compared with the 
asymptotic Pierce theory over a broad range of beam and circuit parameters. A modi- 
fied form of the asymptotic Pierce theory, which did not include the assumption that 
Pierce’s gain parameter was very small, was used. Agreement in this comparison 
was also quite good. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 18, 1977, 

506-20. 
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APPENDIX A 
SYMBOLS 

Pierce launching-loss parameter, defined in eq. (89) 
constant of integration in eq. (22) 

complex amplitude coefficient for induced backward wave 

complex amplitude coefficient for induced forward wave 

unit vector in the z direction 

Pierce gain parameter, defined in eq. (92) 

constant of integration in eq. (22) 

constant of integration in eq, (22) 

Pierce velocity parameter, defined in ref. 12 

Pierce interaction parameter, defined in eq. (90) 

parameter defined in eq. (29) 

parameter defined after eq. (16) 

velocity of light in free space 

Pierce loss parameter, defined in ref. 12 

energy 

electric field due to current source J„ 

a 

amplitude parameter defined in eq. (19) 
axial component of electric field 
electric field 

backward-propagating, source-free electric field 

forward-propagating, source-free electric field 

force parameter defined in eq. (71) 

force parameter defined in eq. (73) 

total axial force 

total axial force on i'^^ disk 

axial component of force associated with k*^” cavity 
axial force on a beam disk due to radio-frequency fields 
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repulsive space- charge force between two beamdii^s 
frequency program input parameter) 
gain (see eq. (88)) 

magnetic field due to current source 
magnetic-field strength 

magnetic- field strength associated with source-free backward wave 
magnetic-field strength associated with source-free forward wave 
beam current (program input p arameter) 
modified Bessel function of order s and argument x 

current density 
current source 

Bessel function of order s and argument x 
axial component of current 

Fourier component of the beam current density 
current 

summation index 

result of surface integration, defined in eq. (34) 

Pierce interaction impedance, defined in ref. 12 
eigenvalue defined by eq. (23) 
periodic length of a cavity 

periodic length of k*^” cavity (program input parameter) 
length of tube 
half of gap length 

half of gap length of k cavity program input parameter) 
electric-field shape parameter used in eq. (16) program input parametei) 
rest mass of a beam disk 
rest mass of a particle 

length of a tube in units of beam wavelengtti 
number of beam radii program input parameter) 
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number of cavities (program input parameter) 

number of disks per beam wavelength program input parameter) 

number of nodes per cavity program input parameter) 

summation index 

a unit vector normal to a surface 

power added to wave by current source 

SL 

fVi 

backward wave power at k cavity 
beam kinetic power at k*'“ cavity 

fVi 

forward wave power at k*^" cavity 

voltage- jump power at k^“ cavity 
th 

power at k cavity 
propagation loss at k^' cavity 

radio-frequency input power (program input parameter) 
summation index 

Pierce space-charge parameter, defined by eq. (93) 
charge of a beam disk 

space- charge reduction factor, defined in ref. 12 
radius 

beam tunnel radius 

beam radius 
th 

i*^" beam radius (program input parameter) 
beam radius in k*' cavity 
a surface 

back surface of a section of a slow-wave structure (see fig. 2) 

front surface of a section of a slow-wave structure (see fig. 2) 

outer, cylindrical, conducting surface of a section of a slow-wave struc- 
ture (see fig. 2) 

summation index 

period of a radio-frequency cycle (l/f) 
transformation defined by eq. (5) 



t 


time 


t“(z) 


u 


U« 


u„ 




V 


'B,k 


B,k 


AV- 

^F,k 

AV 


F,k 


J,k 


Vi, 


AV,. 


V. 


'^0,B 

^0,F 

■if 

sgn(x) 

yn,k 

Z 

Zk 

Z 

Re ('2) 


Zj(t) 


thickness of a beam disk 

arrival time of i “ disk in m*” radio-frequency cycle at axial position z 
velocity 

velocity of i^^ disk 
axial component of velocity 
initial beam velocity 
voltage or potential 
backward voltage on k “ gap 

backward voltage on k^^ gap induced by the beam current in k^^ cavity 
forward voltage on gap 

fVi fVi 

forward voltage on k*^^ gap induced by the beam current in k*"” cavity 

voltage- jump voltage 

voltage-jump voltage on k gap 
til 

voltage on k “ gap 

voltage on k^^ gap induced by the beam current in k*"^ cavity 

beam voltage (program input parameter) 

gap voltage associated with source-free backward wave 

gap voltage associated with source-free forward wave 

volume 

sign of X (positive if x ^ 0; negative if x < 0) 
measure of axial position defined in eq. (69) 
interaction impedance 

interaction impedance of k™ cavity program input parameter) 
axial position 

real part of a complex parameter z 
til 

axial position of disk at time t 

axial position referred to center of k^^ cavity 
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\k 

^n+i,k 

^0 

a 


(«% 

]3 

^3L)k 

0SiLA)k 


y 


n 


°s0 

^0 

M 

Vs 

P,P0 

a 

<P,(pQ 

00 


fh 

length of integration interval for cavity 
position of node in k^^ cavity 

til 

position of n antinode in k “ cavity 
axial separation of two beam disks 
propagation loss per unit length 
power-balance parameter defined by eq. (52) 

propagation loss across k''“ cavity in backward direction (program input 
parameter) 

propagation loss across k*"^ cavity in forward direction (program input 
parameter) 

propagation loss across k*'^ cavity 

phase shift per unit length (axial wave number) 

til 

phase shift per unit length of n space harmonic 
phase shift across k^^ cavity 

parameter proportional to phase shift across k^^ cavity (program input 
parameter) 

radial wave number of n^^ space harmonic 
solution of eq. (94) or (95) 
parameter defined by eq. (11) 

Pierce propagation constants (see eq. (89)) 
permittivity of free space 
beam wavelength 

normalized electric-field shape parameter 

parameter defined by eq. (11) A r is the p*"^ zero of Jg(x)) 

radial coordinates 

charge density 

azimuthal coordinates 

angular frequency (27rf) 

beam plasma frequency (see ref. 12) 
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APPENDIX B 


DERIVATION OF FOURIER COMPONENT OF BEAM CURRENT DENSITY 


In the section on program structure, 3^(z) was defined (eq. (74)) as 


y . 

^ T 4-t 


J(z,t')e"^*^^* dt» 


(Bl) 


We proceed to obtain an ejqiression for Jj^(z) which is convenient for computation. Let 

t™(z) be the time of arrival at z of the disk in the m*'^ radio-frequency cycle. 

We let m = 0 for the radio-frequency cycle considered in the program. To conform 
with previous notation, we will henceforth omit the 0 superscript when referring to 
the cycle considered in the program; that is, t. (z) = t.(z) . Then, for any integer m, 

t“(z) = t.(z) - mT (B2) 

Consider J(z,t) as a function of t for some fixed z (see fig. 9). At t=t?^(z), 

the value of J(z,t) is cru.(z), where a is the charge density of a disk. Furthermore, 

^ in 

J(z, t) has this value for a time Interval centered about t. (z), where the length of the 
interval is tjy'u.(z). We can express this partial contribution to J(z, t) by using the 
unit step function U in the following way: 


au. (z) ( U 



t - t“(z) 




- U«i 


J 2uj(z) 


t - t“(z) 


‘^D 


2u.(z) 


For more concise notation, we define the function h as 


h" 


t - t“(z) 


■'D 


U< 


J 2u.(z) 


J 2u.(z)J |L 


■^D 


2u.(z) 


The partial contribution to J(z, t) can now be written as 


au.(z)h>| t - t^(z) 




J 2u.(z) 


(B3) 
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The total contribution to J(z, t) is obtained by summing over all disks and all cycles. 



By the periodicity of the motion, there will be disks passing the position z in the 
time interval from t-T to t. In general, the N , disks do not all come from the same 
radio- frequency cycle. Let m. be the cycle number of the disk. From equations 
(B2) and (B5) we obtain 


N 


d 


J 


1 





(B6) 


where t.^^ and t -2 


are given by 


til = *i<^) - “iT - 


“D 


2u.(z) 


(B7) 


ti2 = t.(z) - m.T + 


“D 


2u.(z) 


(B8) 


Evaluating the integral in equation (B6) yields 
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i=l 


sm 


at 


cotj 

[2Uj(z)| 


D r~,,t ~i 
2u.(z) 


-io) 

e 


[tj(z)-mjT] 


Since com.T = 27rm., we have 



(B9) 


(BIO) 


The charge density a can be expressed as 


a 


q 


Combining equations (BIO) and (Bll) yields 



sm 


Ci)t- 


D 


2Uj(z) 


tot 


D 


2uj(z) 


-icot.(z) 
e ^ 


(Bll) 


(B12) 
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TABLE 1. - PROPAGATION LOSS PER CAVITY 



TABLE II. - PHASE SHIFT PER CAVITY 




Frequency, f, 

GHz 


■ 

12.05 

12.10 

12.15 

12.20 

Phase shift per cavity, 

(/3^L/ir)j^, rad 

1 

1.48 

1.50 

1.51 

1.53 

1.54 

2 

1.32 

1.34 

1.36 

1.39 

1.41 

3 to 9 

1.16 

1.20 

1.23 

1.27 

1.30 

10 

1.32 

1.34 

1.36 

1.39 

1.41 

11 to 13 

1.48 

1,50 

1.51 

1.53 

1.54 

14 

1.32 

1.34 

1.36 

1.39 

1.41 

15 to 25 

1.16 

1.20 

1.23 

1.27 

1.30 

26 

1.32 

1.34 

1.36 

1.39 

1.41 

27 to 29 

1.48 

1.50 

1.51 

1.53 

1.54 

30 

1.32 

1.34 

; 1.36 

1.39 

1.41 

31 to 39 

1.16 

1.20 

! 1.23 

1.27 

1.30 

40 to 42 

1.16 

1.21 

1 1.25 

1.29 

1.32 

43 to 47 

1.16 

1.20 

1.25 

1.28 

1.31 

48 to 56 

1.17 

1.21 

1.24 

1.27 

1.30 

57 

1.32 

1.34 

1.36 

1.39 

1.41 ; 

58 

' 1.48 

1.50 

1.51 

1.53 

1.54 





TABLE m. - INTERACTION IMPEDANCE 


Cavity 


Frequency, f, 

GHz 



12.00 

12.05 

12.10 

12.15 

12,20 


Interaction impedance, Z, ohms 

1 

400 

400 

400 

400 

400 

2 

850 

800 

780 

750 

730 

3 

1550 

1450 

1375 

1280 

1230 

4 to 9 

1990 

1830 

1725 

1605 

1530 

10 

1550 

1450 

1375 

1280 

1230 

11 

850 

800 

780 

750 

730 

12 and 13 

400 

400 

400 

400 

400 

14 

850 

800 

780 

750 

730 

15 

1550 

1450 

1375 

1280 

1230 

16 to 25 

1990 

1830 

1725 

1605 

1530 

26 

1550 

1450 

1375 

1280 

1230 

27 

850 

800 

780 

750 

730 

28 and 29 

400 

400 

400 

400 

400 

30 

850 

800 

780 

750 

730 

31 

1550 

1450 

1375 

1280 

1230 

32 to 39 

1990 

1830 

1725 

1605 

1530 

40 to 42 

2390 

2015 

1795 

1630 

1535 

43 

2300 

1900 

1710 

1550 

1475 

44 

2150 

1800 

1620 

1470 

1400 

45 

2000 

1700 

1530 

1390 

1330 

46 

1850 

1600 

1440 

1310 

1250 

47 to 51 

1700 

1500 

1350 

1255 

1170 

52 

1600 

1350 

1240 

1150 

1150 

53 

1400 

1200 

1130 

1050 

1060 

54 and 55 

1295 

1115 

1015 

930 

960 

56 

1000 

875 

810 

750 

700 

57 

550 

520 

500 

475 

460 

58 1 

300 

300 

300 

300 

300 


TABLE IV. - CAVITY DIMENSIONS 


Cavity 

Peilodic length, 
cm 

Gap length, 
cm 

1 to 39 

0.3175 

0.0965 

40 to 42 

.3048 

.0914 

43 

.2985 

.0899 

44 

.2921 

.0879 

45 

.2794 

.0843 

46 

.2667 

.0808 

47 to 51 

.2540 

.0772 

52 

.2413 

.0739 

53 

.2286 

.0701 

54 to 58 

.2159 

.0660 


2fK 




: 2r, 


If T TP 


21 , 


Figure 1. - Geometry of coup led- cavity traveling- wave tube model. 




Figure 2. - Geometry of slow- wave structure. 
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Pierce velocity parameter, 
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Figure 3. - Power balance in a cavity. 




(b) Pierce interaction parameter, C, 0.15; Pierce space- 
charge parameter, QC. -0; number of cavities, N(<. 12. 

Figure 4. - Comparison of asymptotic Pierce theory with coupled- 
cavity program. 
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O Asymptotic Pierce theory 

□ Coupled-cavity program 

Pierce velocity parameter 



Beam voltage. Vq. kV 

(d) Pierce interaction parameter. C, 0.05; Pierce space- 
charge parameter. QC. 0.25; number of cavities. N^, 40. 

Figure 4. - Concluded. 
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Figure 5. - Normalized circuit velxity as function of cavity number. 
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12.00 12.05 12.10 12.15 12.20 12.25 12.30 

Frequency, GHz 

Figure 6. - Small-signal gain versus frequency (OdBm drive). 
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Figure 7. - Input power required to saturate tube versus 
frequency. 



Figure 8. - Saturation output power versus frequency. 


J(z.t) 



Figure 9. - Current density J(z. t) as a function of t for some 
fixed z. 
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